PITHA 97/42, HLRZ1997_65 

Strongly coupled U(l) lattice gauge theory 
as a microscopic model of Yukawa theory 



W. Franzki and J. Jersak 
Institut fiir Theoretische Physik E, RWTH Aachen, Germany 
(February 1, 2008) 



> 
o 

o 

(N 

(N 
> 
00 

m 
o 



l> ■ 

Oh, 
^ ■ 



X 



Dynamical chiral symmetry breaking in a strongly coupled 
U(l) lattice gauge model with charged fermions and scalar is 
investigated by numerical simulation. Several composite neu- 
tral states are observed, in particular a massive fermion. In 
the vicinity of the tricritical point of this model we study the 
effective Yukawa coupling between this fermion and the Gold- 
stone boson. The perturbative triviality bound of Yukawa 
models is nearly saturated. The theory is quite similar to 
strongly coupled Yukawa models for sufficiently large coupling 
except the occurrence of an additional state - a gauge ball of 
mass about half the mass of the fermion. 



The question whether the Higgs- Yukawa mechanism of 
symmetry breaking and particle mass generation in the 
contemporary high energy physics might be an effective 
theory to some more fundamental renormalizable quan- 
tum field theory with dynamical symmetry breaking has 
been raised many times. It is of particular interest with 
respect to the large top quark mass, stimulating e. g. the 
idea of topcolor. The pursuit of this question leads fre- 
quently to devising strongly coupled gauge theories be- 
yond the standard model with speculative nonperturba- 
tive dynamics. It might be inspiring to have prototypes of 
such "microscopic" gauge models of Yukawa theory with 
a solid and accessible dynamical framework. Numerical 
lattice methods are justified for this purpose provided a 
lattice model can be found with required basic dynamical 
properties, though neglecting various phenomenological 
aspects. In particular, these properties should be pre- 
served when approaching the continuum limit, i.e. the 
model should be nonperturbatively renormalizable. 

A promising lattice model has been suggested in 
Ref . . This "xt^ (p model" consists of a charged fermion 
field X with strong vectorlike coupling to compact U(l) 
gauge field U, and a scalar field cj) of the same charge. 
Here we summarize the results of our extensive system- 
atic investigations of this model in four dimensions by 
means of numerical simulations with dynamical fermions. 
Though we cannot decide by currently available means 
whether the model is renormalizable, we find some en- 
couraging properties. The detailed account is given in a 
parallel paper |^ and in some previous works . 

The same model has been investigated also in lower 
dimensions. There is little doubt that the model 
defines a renormalizable quantum field theory in two di- 



mensions, belonging to the universality class of the chiral 
Gross- Neveu model [Q . Also first results in three dimen- 
sions suggest renormalizability 

Finding a "promising" lattice model meets two major 
difficulties. Strongly coupled lattice gauge theories with 
fermions exhibit frequently the required dynamical chiral 
symmetry breaking. However, this phenomenon comes 
mostly along with undesirable confinement of fermions 
acquiring mass. Furthermore, at strong coupling, it is 
difficult to find a second order phase transition required 
for an approach to continuum. 

In the x?7(/) model, the scalar field cj) helps to solve 
both problems. Firstly, it shields the fermion charge and 
gives rise to an unconfined, i.e. physical massive fermion 
F = (/)^X in the phase with chiral symmetry broken dy- 
namically by the gauge interaction (Nambu phase). Sec- 
ondly, the scalar suppresses this symmetry breaking and 
at sufficiently strong gauge coupling induces a second or- 
der transition to a chiral symmetric (Higgs) phase, thus 
opening a way to continuum. 

We find that one particular point of the phase transi- 
tion line, the tricritical point E, represents a new kind 
(new universality class) of dynamical symmetry break- 
ing mechanism in strongly coupled gauge theories with 
fermions in four dimensions. We present some evi- 
dence that a strongly coupled effective Yukawa theory 
results. It describes the interaction between the com- 
posite fermion F and some XX "mesons" . At this point 
the model could serve as a microscopic model of Yukawa 
theory. But in addition, the mechanism produces also a 
scalar gauge ball. 

To explain these findings, we now summarize the most 
relevant features of the model (more details can be found 
in 1^]). The model is defined by the action 

S = + Sjj + 5*0 , 



Su = -/3^cos(ep) 



5*0 = -K^^^^{(j)lUx,fj.(f>x+i^L + h.c). 
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FIG. 1. Phase diagram of the xU4i model in the chiral 
limit, mo = 0. The /3 = limit corresponds to the NJL 
model. The NE line is a line of second order phase transi- 
tions. Other lines are lines of first order transitions. In the 
Nambu phase, chiral symmetry is broken dynamically and 
the fermion F = (f)^x is massive. In the Higgs phase, F is 
massless. E is the tricritical point where the scaling behavior 
is different from the NJL model. The Coulomb and Higgs 
phases extend until /3 = oo. The data points are positions of 
the phase transitions extrapolated to the infinite volume. 



Here 6p G [0, 27r) is the plaquette angle, i.e. the argu- 
ment of the product of U(l) gauge field link variables 
Ux,ti along a plaquette P. Taking 8p = o^gF^^, where a 
is the lattice spacing, and /3 = one obtains for weak 

coupling g the continuum gauge action Su = j J d'^xF^^. 
The staggered fermion field x has (real) bare mass aniQ 
in lattice units. It leads to four fermion species in the 
continuum limit. Sign factors r]^^ are standard for stag- 
gered fermions. The complex scalar field is constrained, 
\(f)\ — 1. As its "hopping parameter" k increases, it drives 
the model into the usual Higgs phase. 

We note that there is no Yukawa coupling between x 
and (j), as both fields have the same charge. The model 
has U(l) global chiral symmetry in the limit of vanishing 
fermion bare mass in physical units, rriQ = 0. This is the 
case we are really interested in. The numerical simula- 
tions have to be carried out at nonvanishing mo, however, 
and an extrapolation to the chiral limit performed. For 
amo 0, the chiral transformation relates the model at 
±amQ. 

The phase diagram at toq = is shown in Fig. |l|. One 
limit case is k = 0, corresponding to the massless com- 
pact lattice QED. Another important limit case is /3 = 0. 
Here the model can be rewritten exactly as the lattice 
Nambu- Jona-Lasinio (NJL) model with the critical 
point N. The strong four-fermion coupling of that model 
corresponds to small k. The Nambu phase thus connects 
the broken chiral symmetry phases of both these limit 
models. Only a part of its boundary, the NE line, repre- 
sents second order phase transitions. 



Point E is far away from any limit case and it does not 
appear to be accessible by any reliable analytic method, 
neither on the lattice nor in continuum. It is "tricritical" 
because in the full parameter space (including amo) there 
are, apart from NE, two further second order "wing" lines 
entering E from the positive and negative amp directions. 
The existence of a common point E of these three sec- 
ond order lines is neither predicted nor understood. The 
evidence is purely numerical, but quite strong Its 
position is /3e = 0.62(3), ke = 0.32(2). 

The importance of the point E roots in the experience 
from statistical mechanics that a tricritical point belongs 
to a universality class different from that of any of the 
second order lines entering into it |10|. (Basic properties 
of tricritical points relevant in our context are summa- 
rized e.g. in ||ll[].) The usual universality of second order 
lines suggests, and this is supported by our earlier inves- 
tigations that the whole NE line except the point E 
corresponds to the same continuum model as the point 
N, the NJL model. The gauge field is presumably auxil- 
iary and the model is therefore of limited interest there. 
However, the point E is expected to be different, gauge 
field playing an important role. 

To verify this expectation, we have investigated criti- 
cal exponents and spectrum of the model in the vicinity 
of E. The study of exponents uses advanced techniques 
of statistical physics and is described in j^]. Here we 
only mention the found value vt — ^/S of the correlation 
length tricritical exponent. It is different from ~ 1/2 
along the NE line confirming the particularity of E. It 
also differs from the prediction of the classical theory of 
tricritical points usually believed to hold in four dimen- 
sions and predicting = 1/2. This indicates that 
the point E is a tricritical point with important role of 
quantum fluctuations. 

Some insight into the continuum physics, which might 
be obtained at the point E, is provided by the spectrum 
and its scaling behavior. The masses amq are in lat- 
tice units and thus only their ratios are physical. They 
are determined in the Nambu phase without any gauge 
fixing from the correlation functions of various gauge in- 
variant composite operators Q. In this sense the massive 
physical fermion F, as well as other physical states, are 
composite. The interaction between them is due to the 
Van der Waals remnant of the fundamental interactions. 

The fermion mass amp > decreases when the NE 
line is approached from the Nambu phase and there it 
is rather insensitive to the lattice size. It is small in the 
Higgs phase on finite lattices, decreasing with increasing 
lattice size. The data are consistent with expected van- 
ishing of amp in the Higgs phase in the infinite volume 
limit. 

We find several x, x bound states and borrow their 
names from QCD. One of them is the obligatory pseu- 
doscalar Goldstone boson tt with the dependence on amo 
as required by current algebra. The 7r-meson is mass- 
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less in the chiral limit in the whole Nambu phase. It 
would get "eaten" if the global chiral symmetry of the 
XU4> model were gauged. This process would be identi- 
cal with the standard Higgs mechanism. Thus we do not 
need to discuss it and mention it only for completeness. 

We obtain some results for the pion decay constant. 
The data suggest a large ratio j-Kimp increasing when 
the point E is approached. However, the value is sensitive 
to the lattice volume and we cannot yet extrapolate it 
to the infinite volume and continuum limit. If the ratio 
diverges, it would indicate the triviality of the model [ p2| . 
Its current value, fT^/mp ~ 1/3 for amp ~ 0.4 can be 
considered as a lower bound. 

Also the a meson in the antifermion-fermion channel 
is observed. The a mass is quite dependent both on the 
lattice size and bare fermion mass, preventing its predic- 
tion in the continuum limit with mo = 0. Strong lattice 
size dependence of and a mass has been observed and 
explained by means of the Schwinger-Dyson equations in 
the limit case /? = 0, the NJL model, and holds appar- 
ently also in the vicinity of E. 

The p mass is insensitive to the lattice size and scales 
like the fermion mass, with the approximate value mp ~ 
2m p . We cannot distinguish between a bound state and 
a resonance. 

Also a neutral scalar (S'-boson) is seen appearing as 
a composite of 0^ and (p and as a state of pure gauge 
field. In the Higgs phase it corresponds to the Higgs bo- 
son associated with the perturbative Higgs mechanism 
occuring in that phase at large (3. In the Nambu phase 
it is more natural to interpret the S'-boson as a scalar 
gauge ball. The transition between both interpretations 
is smooth, however, and the corresponding channels ap- 
preciably mix in the vicinity of E. 

The mass ams is nonzero in all phases and goes to 
zero when the two wing critical lines are approached for 
any amp. For amo = 0, the mass ams goes to zero 
only at the point E. This illustrates the particular char- 
acter of the tricritical point E. Whereas only one of the 
masses amp or ams vanishes on each of the second or- 
der lines, they both vanish at the tricritical point. As 
their ratio remains finite, both corresponding states are 
present in the continuum limit taken at this point. This 
is the best evidence that the physical content (universal- 
ity class) of the point E differs from that of any of the 
adjacent second-order lines and thus does not correspond 
to the NJL model obtained at N. 

As ams is only moderately dependent on the lattice 
size, we can estimate its scaling behavior when point E is 
approached. This is shown in Fig. ^for k = 0.30, which is 
approximately the K-coordinate of E. The ratio ms/mp 
remains constant when amp decreases, as long as finite 
size effects are negligible. The sudden rise of the data 
at small amp is due to the strong finite size dependence 
of amp in the Higgs phase and shifts correspondingly to 
smaller amp when the lattice volume increases. These 



1,6 
1.4 
1,2 
1,0 

o,a 

0,6 
0,4 
0.2 



/< = 0.30 
am„ = 0.02 



0:6n6 
•:8'24 



0.2 0.4 0.6 0.8 



FIG. 2. Mass ratio ms/niF as function of the fermion mass 
when point E is approached. The data have been obtained 
on 6^16 (empty circles) and 8*^24 (full circles) lattices. 



observations suggest that the value ms/mp ~ 1/2 would 
be obtained in the continuum limit at E. 

In the rest of this letter we concentrate on the renor- 
malized Yukawa coupling ?/r between F and tt. It is ob- 
tained from the three-point function of the correspond- 
ing composite operators and thus can be interpreted as 
an effective Yukawa coupling, which would describe the 
interaction between F and tt in the regime where their 
composite structure can be neglected. Our measurement, 
performed in a similar way as in ||l3|] , is described in detail 
in [^j6| . In spite of the complexity of the corresponding 
expressions, j/r is measurable with good precision. The 
reason is the strict locality of the used operators. An 
attempt to measure also an analogous coupling of the 
S'-boson failed because S is described by extended oper- 
ators. 

The renormalized Yukawa coupling j/r is presented in 
Fig. H (full symbols) for three different amo at k = 0.30 
in the Nambu phase close to E. We show it in dependence 
on amp and compare it with the tree level relation 



(tree) _ amp 
" (XX) 



(1) 



The agreement is so good that the open symbols repre- 
senting (^ in fig. 1^ are nearly invisible. 

The values for j/r are not significantly dependent on 
amo. Therefore we expect that the value for in the 
chiral limit is rather well represented by our results at 
amo = 0.01. In Fig. ^ we show these results for substan- 
tially different gauge couplings. The results close to the 
point E, for j3 — 0.64, n — 0.30, and (3 — 0.55, are very 
similar. The values of 2/r at /3 = (NJL model) are signif- 
icantly below the other data, however. This implies that 
the Yukawa coupling gets stronger as one moves from the 
NJL model to the vicinity of E by increasing (3 at fixed 
amp. 

In the interval j/r ~ 2 — 5, where the Yukawa coupling 
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For j/R < 2, the data are significantly influenced by the 
finite size effects, and we cannot make any conclusion. It 
might be that some deviation from the Yukawa theory 
occurs there. The fact that j/r is stronger in the vicin- 
ity of E than in the NJL model suggests the question 
whether it vanishes in the continuum limit, amp — 0, 
taken there. Is it zero as in Yukawa theory, or could 
triviality be avoided at the tricritical point? A reliable 
extrapolation of ?/r to am i? = at the tricritical point of 
the xU(t> model is a challenge for future investigations. 
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FIG. 3. Effective Yukawa coupling of the n meson to F 
(full symbols), and the tree level relation (|l|) (open symbols) 
as function of amp for different amo at k = 0.30 on 8"^24 
lattice. 
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FIG. 4. yn, of the tt meson for different couplings along the 
line NE. For f3 — 0.55 the tree level definition was used. Most 
data have been taken on a 8^24 lattice (at /? = on 6^16) 
with amo ~ 0.01 (at P = 0.55 is amo ~ 0.02). Shown is also 
the curve of maximal as it results from the perturbative 
expansion first order in a coresponding Yukawa model. This 
curve goes to for amp 0, as expected from triviality. 



is determined reliably, we compare our result with the 
curve of maximal renormalized coupling, resulting from 
the first order perturbative calculation in the Yukawa 
model (triviality bound): 
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yRmax ^ 



^2f3o\\n{amF)\ ' 



(2) 



with Np — 4, the number of fermions including lattice 
doublers. In the lattice Yukawa models this upper bound 
is nearly saturated by the values of obtained at the 
maximal possible bare Yukawa coupling |l^Jl^. As seen 
in fig. ^ also in our model close to E, the data are only 
slightly below the curve (|^) . Thus in this interval, yn has 
the same form as that in the strongly coupled Yukawa 
models. 
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